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OPTICAL SPECTROPOLARIMETRY OF THE GRB 020813 AFTERGLOW 
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ABSTRACT 

The optical afterglow of gamma-ray burst 0208 1 3 was observed for 3 hours with the LRIS spectropolarimeter at 
the Keck-1 telescope, beginning 4.7 hours after the burst was detected by HETE-2. The spectrum reveals numerous 
metal absorption lines that we identify with two systems at z = 1.223 and z = 1.255. We also detect an [O II] 
A3727 emission line at z = 1.255 and we identify this galaxy as the likely host of the GRB. After a correction 
for Galactic interstellar polarization, the optical afterglow has a linear polarization of 1.8-2.4% during 4.7-7.9 
hours after the burst. A measurement of p = 0.80% ± 0.16% on the following night by Covino et al. demonstrates 
significant polarization variability over the next 14 hours. The lack of strong variability in the position angle of 
linear polarization indicates that the magnetic field in the jet is likely to be globally ordered rather than composed 
of a number of randomly oriented cells. Within the framework of afterglow models with collimated flows, the 
relatively low observed polarization suggests that the magnetic field components perpendicular and parallel to the 
shock front are only different by about 20%. 

Subject headings: gamma rays: bursts — polarization 
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1. INTRODUCTION 

It is now generally assumed that the optical afterglows of 
gamma-ray bursts (GRBs) arise from synchrotron emission in 
the post-shocked gas of a relativistic blast wave (see, e.g., Piran 
2000; Meszaros 2002). There is credible evidence that GRBs 
are strongly collimated or beamed explosions (e.g., Rhoads 
1999; Sari, Piran, & Halpern 1999; Frail et al. 2001). Syn- 
chrotron radiation is intrinsically polarized, and the synchrotron 
emission from a jet will have a net polarization if its axis does 
not coincide with the line of sight. A natural consequence 
of models for beamed afterglows is time-variable polarization, 
which will depend on the viewing geometry as well as the de- 
gree of orientation of the magnetic field (e.g., Medvedev & 
Loeb 1999; Gruzinov 1999), up to a maximum polarization of 
- 20% (Ghiseflini & Lazzati 1999; Sari 1999). Optical imag- 
ing polarimetry has been performed for a handful of GRB after- 
glows to search for this signature. The results have been either 
upper limits (Hjorth et al. 1999; Covino et al. 2002a) or detec- 
tions at p « 1 - 3% (Covino et al. 1999; Wijers et al. 1999; Rol 
et al. 2000; Covino et al. 2002b), except for an observation of 
GRB 020405 by Bersier et al. (2002) that detected 9.9% ±1.3% 
polarization 1.3 days after the burst. 

GRB 020813, a typical long-duration burst with a series of 
pulses, was detected by HETE-2 at 2:44 UT on 2002 August 
13 (Villaseiior et al. 2002), and an associated optical transient 
(OT) was promptly identified (Fox, Blake, & Price 2002). The 
rapid localization allowed us to begin spectropolarimetric ob- 
servations less than 5 hours after the burst was initially detected. 
These are the first spectropolarimetric observations obtained for 
a GRB afterglow. Preliminary descriptions of this dataset were 
given by Price et al. (2002) and Baith et al. (2002). 



2. OBSERVATIONS AND REDUCTIONS 

The observations were obtained with the dual-beam LRIS 
spectropolarimeter (LRISp; Oke et al. 1995; Goodrich, Cohen, 
& Putney 1995) on the Keck-I telescope, and a slit width of 1". 
The addition of a new ultraviolet polarizing calibration filter to 
LRISp allowed us for the first time to take advantage of the full 
available wavelength range of LRISp using both the blue and 
red cameras. We used a 400 lines mm"' grism on the blue side 
of the spectrograph and a 400 lines mm"' grating on the red 
side, giving a wavelength scale of 1.09 A pixel"' over 3200- 
5800 A in the blue and 1.86 A pixel"' over 5600-9400 A in 
the red. The red and blue beams were separated with a dichroic 
beamsplitter having 50% transmission at 5700 A. 

We began exposures of the OT associated with GRB 020813 
on 2002 August 13 at 7:23 UT, when its brightness was V w 19 
mag. The total exposure time was 3 hours, broken into individ- 
ual 15-minute exposures. Each pair of consecutive exposures 
gives one of the Stokes parameters q (with the half-wave plate 
oriented at 0° and 45°) or u (22.°5 and 67.°5), and tiie 4 expo- 
sures taken together are used to determine the degree and po- 
sition angle of linear polarization (p and 9). We refer to these 
three sets of four exposures each as the "hour 1", "hour 2", 
and "hour 3" data. The airmass was 1.3-1.5 during the obser- 
vations, and before starting each hour's exposures the spectro- 
graph was rotated to place the slit at the parallactic angle for 
the midpoint of the hour Seeing was ^ Q'.'l, and there were 
intermittent thin clouds during the observations. 

As a probe of the Galactic interstellar polarization (ISP), 
we also observed the A2V star HD 187330, which is located 
64 .'7 from the line of sight to the GRB. Its apparent magnitude 
(V = 9.22 mag) implies a distance of ^ 350 pc; this is suffi- 
ciently distant to probe the majority of the Galactic ISP for a 
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Galactic latitude of -21° (Tran 1995). Distant A-type stars are 
useful as probes of ISP because of their brightness and low in- 
trinsic polarization, although ideally a smaller angular separa- 
tion from the OT would be desirable. 

The data were reduced using standard techniques as de- 
scribed by Miller, Robinson, & Goodrich (1988). The zeropoint 
of 6 was calibrated with the polarized standard star BD+59°389 
(Schmidt et al. 1992). For the unpolarized standard star HD 
14069 we found p<0.1% across the entire optical spectrum. 

3. CONTINUUM AND ABSORPTION LINES 

Figure 1 displays the total flux spectrum of the OT from all 
exposures combined. Numerous metal absorption lines are su- 
perposed on a featureless continuum, as is typical for GRB af- 
terglow spectra (e.g., Metzger et al. 1997; Kulkami et al. 1999; 
Vreeswijk et al. 2001; Jha et al. 2001). Two distinct absorption 
systems are present in the spectrum (Table 2). The individ- 
ual absorption lines in the higher-redshift system have a mean 
z = 1.255. The second system has weaker absorption lines and a 
mean redshift of 1.223. Both systems show lines of Si II A1527, 
C IV A A 1 548 , 1 55 1 , Al III A 1 67 1 , Mg II A A2796 ,2804, and Mg I 
A2853, as well as several Fe II lines including AA2587,2600. 
Additional features from the z = 1.255 system include lines of 
Al III, Mn II, Mg I, Zn II, and Cr II. One emission line is also 
present in the spectrum, centered at 8404.0 A and consistent 
with [O II] A3727 at z = 1 .255. We conclude that the z = 1 .255 
system is the Ukely host of the GRB, since it is the highest- 
redshift absorption system detected and its [O n] emission line 
is an indicator of star-forming activity. 

Ionization of the surrounding interstellar medium (ISM) by 
the GRB could cause temporal variations in the equivalent 
widths of interstellar absorption lines within the host galaxy 
(Pema & Loeb 1998), and we have measured the equivalent 
widths of the absorption lines in the hour 1, 2, and 3 spectra sep- 
arately to search for any variations. The differences in equiva- 
lent widths are consistent with the level of measurement uncer- 
tainty, and we conclude that no variability is detected. Searches 
for ultraviolet absorption-line variability in other GRBs have so 
far yielded nuU results as well (Vreeswijk et al. 2001; Mirabal 
et al. 2002). 

The continuum has a power-law slope of /3 = -1.06 ± 0.01 

(for /,y oc i^^) on the red side. On the blue side, we find a steeper 
slope of /3 = -1 .25 ± 0.03. This apparent steepening toward the 
blue is affected by atmospheric dispersion and is sensitive to the 
centering of the OT within the slit, as well as reddening, so the 
red-side value is a more reliable measurement of the intrinsic 
slope. Levan et al. (2002) reported a shallower optical slope of 
P = -0.8 approximately 98 hours after the burst. 

4. POLARIZATION 

Polarization was measured over five wavelength bins: 3500- 
4500, 4500-5500, 5800-6800, 6800-7800, and 7800-8800 
A. AU measurements were performed using the q and u spectra 
and their associated error spectra, with the results converted to 
p and 9 as the final step. We observed the same three polarized 
standard stars (BDh-59°389, HD 19820, and HD 155528) dur- 
ing both nights of our observing run (August 12 and 13 UT). For 
all three stars and five wavelength bins, the RMS difference in 
the polarization measured from the two nights' observations is 
only 0.04%. The polarized standard stars also have essentially 

^ Covino et al. (2002c) use an ISP correction derived from stars in the same 
reasonable agreement with the polarization we measure for HD 187330. 



constant 6' as a function of wavelength (to within ±1°), and 
agreement to within 1° between the two nights of observations. 
This demonstrates the exceUent stability of the instrument. 

The raw measurements for the OT yield /? = 2.3 — 3.1% with 
9 = 153° - 162°, but a portion of this polarization is caused 
by transmission through regions of aligned dust grains in the 
Galactic ISM. The line of sight to GRB 020813 has a Galac- 
tic reddening of E{B-V) = 0.111 mag (Schlegel, Finkbeiner, 
& Davis 1998). Along most Galactic lines of sight, the max- 
imum interstellar polarization observed is 9% per magnitude 
of E(B-V) (Serkowski, Mathewson, & Ford 1975). The ISP 
probe star has p = 0.67% ± 0.01% at 6» = 167.°0 ± 0.°3 over 
5800-6800 A, falling comfortably within the expected range, 
and its p{X) spectrum is consistent with a Serkowski interstellar 
polarization curve. To correct the polarization measurements of 
the OT, we fitted a low-order polynomial to the q and u spectra 
of the ISP star and subtracted the result from the Stokes param- 
eters of the OT. The ISP-corrected measurements for the OT are 
listed in Table 1; random uncertainties due to photon statistics 
and detector read noise are at the level of 0.1% in p and 1-2° 
in 6. The standard star observations suggest that any additional 
systematic uncertainties should not be time-variable, and are 
likely to be < 0.1% in p and < 1° in 0. 

The ISP-corrected polarization longward of 5600 A is nearly 
constant during the 3 hours at 1 .9-2.2%. On the blue side, how- 
ever, the polarization drops from 2.4% ±0.1% during hour 1 to 
1.8% ± 0.1% in hour 3. This is formally a 4cr detection of po- 
larization variability, but it is somewhat surprising to find sig- 
nificant variations in p only on the blue side. As a check on this 
measurement, the reductions were repeated with different ex- 
traction and background windows, and repeated again using en- 
tirely different reduction software; in all cases the results were 
consistent with the original measurements to within 0.1% in p. 
We conclude that the variability in p on the blue side is likely 
to be genuine, although we caution against over-interpretation 
of this result since the measurement is not repeatable. 

Imaging polarimetry obtained on the following night at the 
VLT by Covino et al. (2002c) demonstrates a much more dra- 
matic level of variability in p over the next 14 hours. Covino 
et al. find a raw polarization of p = 1.17% ±0.16% at 6* = 
158° ± 4°, and ISP-corrected values ofp = 0.80% ± 0.16% at 
9 = 144° ±6°.^ If we apply our ISP correction to their raw 
measurement we obtain p = 0.57% at ^ = 147°. The change in 
p between our hour 3 data and the VLT observation is signif- 
icant at the 7(7 level. This large drop in p demonstrates that 
most of the polarization we detect in the Keck data must be in- 
trinsic to the OT itself, and not merely the result of transmission 
through dust within the host galaxy or the intervening galaxy at 
z= 1.223. 

The ISP-corrected polarization detected by Covino et al. 
(2002c) is sufficiently low that in principle it could arise from 
the interstellar medium in the host galaxy or perhaps even in 
the intervening z = 1.223 system, although dust grains have a 
low polarization efficiency in the ultraviolet (Martin, Clayton, 
& Wolff 1999). Furfliermore, if GRBs result from the deafli of 
massive stars, one can expect many GRBs to occur in dusty and 
gas-rich star-forming regions, and in extreme cases Type II su- 
pernovae can have optical ISP of several percent from the host 
galaxy (Leonard et al. 2002). However, Chandra observations 
have indicated a low absorbing column of A^h = 7 x 10^'' cm"^ 

ig field as the OT. Then- V-band ISP correction is p = 0.59% a<L0 = 178°, in 
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that is consistent with the expected Galactic absorption along 
the line of sight (Vanderspek et al. 2002). Also, the lack of 
any large change in accompanying the variation in p suggests 
an intrinsic origin for the majority of the detected polarization. 
Detection of a continued late-time decline in p would provide a 
definitive test. 

5. DISCUSSION 

Our observations, combined with those of Covino et al. 
(2002c), clearly show variability in the degree of polarization, 
with a nearly constant position angle. This already places 
strong constraints on the polarization models. Random effects, 
such as emission from a small number of randomly oriented 
magnetic field cells (Gruzinov & Waxman 1999) or ampUfica- 
tion of the emission from one such cell by microlensing (Loeb 
& Perna 1998), would predict that p and 6 should change on 
the same timescale. 

In principle, the jet geometry, together with an assumption of 
some level of orientation of the amplified magnetic field with 
respect to the shock, could predict the degree of polarization as 
function of time. Models based on this scenario (Ghisellini & 
Lazzati 1999; Sari 1999) show that the maximum degree of po- 
larization occurs around the time of the jet break (fjet), when the 
opening angle of the jet is equal to the inverse Lorentz factor. 
However, at that time the emission is observable from a signif- 
icant fraction of the jet. This introduces order-unity uncertain- 
ties in the polarization predictions since the energy distribution 
profile as function of distance from the axis of the jet is un- 
known. In addition, a detailed description of the hydrodynamic 
evolution is not yet available; only the asymptotic scalings, well 
before and well after fjet, are understood. Thus, it is not possi- 
ble to make detailed, quantitative comparisons of polarization 
models with observations. 

Despite these uncertainties, we can make general compar- 
isons of the observations with the model proposed by Sari 
(1999), which predicts sharp features in the polarization light 
curves, although more realistic assumptions on the hydrody- 
namics of jet spreading, and on the energy distribution within 
the initial jet, are expected to produce smoother features than 
those predicted by this model. This model has three parame- 
ters. The line-of-sight offset to the jet determines the shape of 
the curve; t^^t sets the overall timescale; and an overall scale 
factor for p is set by a function of the ratio of the magnetic field 
components parallel to and perpendicular to the shock front, 
which is assumed to be constant. We denote the offset in units 
of the opening angle of the jet by < < 1 . (For (j)>\ the jet 
is not visible initially, corresponding to an "orphan" afterglow.) 
Small offsets give rise to a single polarization peak, while large 
offsets yield three peaks, where the middle one has position an- 
gle perpendicular to the other two. More realistic jet models 
may smooth the curve and eliminate this middle peak. We note 
that the model of Ghisellini & Lazzati (1999) predicts only two 
peaks since their calculations did not include the spreading of 
the jet. 

Figure 2 demonstrates that the jet model parameters can be 
chosen so as to find at least a qualitative agreement with the 
data. In the figure, our measurements for the 5800-6800 A 
polarization, and the V-band measurement from Covino et al. 
(2002c), are compared with two sample models. Initial analysis 



of the optical hght curve suggested fjet « 3.5-5 hours (Bloom, 
Fox, & Hunt 2002; Gladders & Hall 2002), although subse- 
quent analysis (to be presented elsewhere) indicates that the 
light curves may be consistent with a later value of ?jet w 10 - 1 1 
hours. Given this uncertainty in deriving fjet from the pho- 
tometric data, our approach is to determine whether there are 
any values of ijet that produce polarization light curves consis- 
tent with the observations. The first model has a small offset 
4) = 0.3 and predicts one peak of polarization. (Models with 
< ^ < 0.35 would give similar results.) This requires fjet to 
occur during the Keck observations, at ~ 6 hours. However, 
this model somewhat underpredicts the decay of polarization 
between our measurements and those of Covino et al. A bet- 
ter fit is obtained for a medium-offset model with cf) = 0.55 and 
tjet « 16 hours, where the polarization drops to zero, or perhaps 
even rotates the position angle by 90°, between the Keck and 
VLT observations. Under both interpretations, the overall level 
of polarization is less than the maximal possible level of 20% . 
This implies that the magnetic field perpendicular to the shock 
differs only by 15-20% from that parallel to the shock front. 
Models with other values of have difficulty matching both the 
near-constancy of p during the Keck observations (for A > 5800 
A) and the large change in p detected by Covino et al. (2002c). 

The polarization behavior of GRB 020813 appears rather 
similar to that of GRB 990712, which varied between p = 3.0% 
at 0.4 days after the burst and 1 .2% at f = 0.8 days, at nearly con- 
stant position angle (Rol et al. 2000). We note that our models 
with (j) = 0.3 or 0.55 could also be nearly consistent with the 
GRB 990712 data within its uncertainties. Still, as noted by 
Rol et al., the lack of clear evidence for a jet break for that 
event casts some doubt on this interpretation. 

Finally, we mention that the hints of wavelength dependence 
to the polarization of GRB 020813, which appear to be time- 
variable, can not be accommodated by these models. A spectral 
break passing through the optical band at the time of our three 
measurements could produce such a signature, but there is no 
evidence for such a break. 

To date, polarimetric observations of GRB afterglows have 
had rather sparse temporal coverage. Further progress will re- 
quire continuous polarization monitoring of bright afterglows, 
initiated as early as possible and extending well beyond fjet- 
Also, our data highUght the need for color information in future 
polarization measurements. Obtaining such observations would 
pose some logistical difficulties, but the potential for improved 
understanding of the physical structure of the synchrotron jets 
from GRBs would make this effort extremely worthwhile. 
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Table 1 



POLARIZATION MEASUREMENTS 



Hour 


Wavelength 
(A) 




P 

(%) 






n 


1 


3500-4500 


2, 


.4±0. 


1 


161 ±1 




4500-5500 


2. 


.4±0. 


,1 


159 ±1 




5800-6800 


2, 


,1±0. 


1 


158 ±1 




6800-7800 


2, 


,1 ±0. 


,1 


155 ±1 




7800-8800 


2, 


,1±0. 


1 


155 ±1 


2 


3500^500 


1. 


.8±0. 


1 


160 ±1 




4500-5500 


2, 


,0±0. 


1 


155 ±1 




5800-6800 


2. 


,0±0. 


1 


151±1 




6800-7800 


2, 


,0±0. 


,1 


150 ±1 




7800-8800 


2, 


,1±0. 


1 


151±2 


3 


3500-4500 


1, 


,8±0. 


1 


153 ±1 




4500-5500 


1, 


,8±0. 


,1 


149 ±1 




5800-6800 


2, 


,1±0. 


,1 


156 ±1 




6800-7800 


2, 


,2±0. 


,1 


153 ±1 




7800-8800 


1, 


.9±0. 


1 


149 ±1 



Note. — The results given above have been cor- 
rected for ISP using observations of the A2V star HD 
187330, as described in the text. The midpoints of the 
hour 1, 2, and 3 observations were at 7:54, 9:01, and 
10:06 UT, respectively. 
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Table 2 

Absorption Lines in the GRB 020813 Spectrum 



Ahelio 


Ion 


Arest 


z 


Wx (rest) 


(A) 




(A) 




(A) 




TT 


1 jZO. / 


1 


U.40 ± U. IZ 


J44Z.J 


L. IV Diena 


1j4o.Z, IjjU.o 


1 ooo 
l.ZZZ 


Z.Uj ± U.UO 




1 Ci TT 
+ Ol 11 


1 jZO. / 


1 o<< 
l.ZjJ 












i . / i it U. lo 


'XA CO /I 


TA7 

IV 


1 j4o.Z 




i .J J ± U.UJ 


Q/IQ^ 1 


L. IV 


1 jjU.o 


l.Zj4 


l.jl itU.UO 


JOZO.J 


Tja TT 

re 11 


lOUo.J 


1 .Zj4 


l.Uj ±U.U/ 


3 / 14.0 


A 1 TT 
Al 11 


10 /U.o 


i .ZZj 


U.oo ± U.U4 


J /oo.y 


A 1 TT 
Al 11 


10/U.o 


1 

l.ZjJ 


1.4o ± U.U4 


A(YT1 ^ 
4U / / .J 


Ci TT 

ol 11 


1 ens n 

loUo.U 


1 

l.ZjJ 


n T 1 4- n no 
u. / 1 ± u.uz 


/1 1 A 
41oZ.O 


A 1 TTT 
Al 111 


1oJ4. / 


1 .ZD J 


yj.y L ± u.uz 


4ZUU.y 


A 1 TTT 
Al 111 


loOZ.o 


i.ZjJ 


U.JO ± U.Uj 


4joy.u 


ivig 1, An 11 , d 11 


ZUZO.J, ZUZO.l, ZUZO.J 


i .Zj J 


U.Oi ± U.Uj 


40J0.J 


TT 

KS 11 


ZUjO.j 


1 

l.ZjJ 


U.jU ± U.Uj 


40JU.0 


f^f TT 7-r\ TT 

L^r 11, Zin 11 


zuoz.z, zuoz. / 


1 

l.Zjj 


U.44 ± U.UJ 


40 J O.J 


L,r 11 


zuoo.z 


1 OSS 


U. 1 / it U.Uj 


4oo4.Z 








U.JO ± U.U/ 




T7a TT 

re 11 


zz4y.y 


i .zzz 


U.jj ± U.Uj 


jU/Z.j 
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Note. — Observed wavelengths are heliocentric, and converted to vacuum. 
Rest wavelengths are given as vacuum wavelengths. Absorption-line equivalent 
widths (Wx) and uncertainties were measured with Gaussian fits using the IRAF 
SPLOT task, and have been transformed to the rest frame. 

''For unidentified absorption Unes, equivalent widths are given in the observed 
frame. 
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Fig. 1. — Total flux spectrum of GRB 020813, from all LRISp exposures combined. Absorption features from the z = 1.255 system are labelled above the 
spectrum, and features from the z = 1.223 system are labelled below. Residuals from the removal of the telluric A and B bands are marked with ©. 
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Fig. 2. — Polarization behavior of GRB 020813 compared with model calculations. Squares denote the Keck observations in the 5800-6800 A band, and the 
open circle is the VLT measurement from Covino et al. (2002c). The solid curve is the model with (j> = 0.3 and tj^t = 6 hours, and the dashed curve is the model with 
4> = 0.55 and tj^t = 16 hours. Negative polarizations in this plot indicate a rotation of 90° in the position angle relative to positive polarizations. For consistency, we 
have applied our ISP correction to the VLT data point. 



